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Abstract 

This paper is devoted to provide a theoretical underpinning for 
ensemble forecasting with rapid fluctuations in body forcing and in 
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boundary conditions. Ensemble averaging principles are proved under 
suitable 'mixing' conditions on random boundary conditions and on 
random body forcing. The ensemble averaged model is a nonlinear 
stochastic partial differential equation, with the deviation process (i.e., 
the approximation error process) quantified as the solution of a linear 
stochastic partial differential equation. 

Short Title: Ensemble Averaging under Random Boundary Con- 
ditions 

Key Words: Multiscale modeling, ensemble averaging, random 
partial differential equations; stochastic partial differential equations; 
random boundary conditions; martingale. 

Mathematics Subject Classifications (2010): 60Hf5, 60F10, 
60G17 

1 Motivation 

A complex system often involves with multiple scales, uncertain parameters 
or coefficients, and fluctuating interactions with its environment. Ensemble 
forecasting for such a complex system is a prediction method to obtain col- 
lective or ensembled view of its dynamical evolution, by generating multiple 
numerical predictions using different but plausible realizations of a model 
for the system. The multiple simulations are generated to account for errors 
introduced by sensitive dependence on the initial/boundary conditions and 
errors introduced due to imperfections in the model [3 [8]. 

In order to better understand the theoretical foundation of ensemble fore- 
casting, we consider a system modeled by partial differential equations (PDEs) 
with fast oscillating random forcing in the physical medium or on the phys- 
ical boundary, and show that ensemble averaged dynamics converges to the 
original dynamics, as a scale parameter tends to zero. 

Some relevant recent works [H [3], [Qj [T3] are about averaging or homoge- 
nization for random partial differential equations (random PDEs) with fast 
oscillating coefficients in time or space. Different from the method in the 
above mentioned works, we present a more direct approach, in order to derive 
an averaging principle and deviation estimates for PDEs with random oscil- 
lating coefficients and random oscillating boundary conditions. We previ- 
ously studied [12] stochastic partial differential equations (stochastic PDEs) 
with perturbed white noise dynamical boundary conditions which are mea- 
sured by a small scale parameter e > 0. In that case, the effectively reduced 
model does not capture the influence of the random force on boundary. 
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In the present paper, we consider the following PDE with a random os- 
cillating body forcing and/or a small random oscillating boundary condition, 
for a unknown random field u € (x, t, u) 



Here x G (0,/), / > 0, t > 0, e is a small positive scale parameter, and u is 
in a sample space Q. A probability P with a a— algebra J 7 is defined on this 
sample space. The mathematical expectation with respect to P is denoted 
by E. We often suppress the u— dependence for notational clarity 

We prove an ensemble averaging theorem (Theorem [T] in §2j) for the ran- 
dom system (JXj) — (Ej) , i-e., we obtain a stochastic model for u e as e — > 0. It 
turns out that the random boundary condition appears as a white noise on 
the dynamical field equation for u, as e — > . The ensemble averaged model 
is a stochastic partial differential equation (stochastic PDE) for u, instead of 
a random PDE, with a homogenous boundary condition. When the random 
boundary condition is absent, we further show that the deviation process 
(i.e., approximation error process), u £ — u, can be quantified as the solution 
of a linear stochastic PDE (Theorem [2] in §3]). 

On the technical side, in order to pass the limit e — > , we first prove 
the tightness of the distribution of {u e } , so we just consider (u e , cp) for every 
bounded continuous function tp with compact support. Then, in §2.21 we 
construct a process Ail, which is a martingale by Ethier and Kurtz's result jU 
Proposition 2.7.6]. This construction is very direct [3]. By passing the limit 
e — > in Ail, we obtain the stochastic PDE satisfied by the limit u of u e . 
This method is also applied to show that the deviation process, u e — u, is the 
solution of a linear stochastic PDE; see § [3j 

Note that we take A/e as the intensity scale for the noise boundary con- 
dition. This is for simplicity. In fact our approach can also treat the case 
€2 A1 f(t/e a ) , with < a < 2. A similar case is also discussed in [3]. But 
the case a > 2 is more singular, one should consider the limit of e^" 1 ^ as 



This paper is organized as follows. After recalling some basic background, 
we prove an ensemble averaging theorem for a random PDE system with a 
random boundary condition and with a random body forcing, in §0 and 
further characterize the deviation process in §|3j 



u t = u lx + g(t/t,u e ,uj) , u e (x,0) = u (x), 
u e (0,t) = V£/(*/e,w), u*(l,t) = 0. 



(1) 
(2) 



e ->• 0. 
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2 Ensemble averaging under small fast oscil- 
lating random boundary conditions 

We consider the random PDE system ([I])-([2]). Consider the Hilbert space 
H = L 2 (0,l) with the usual norm || ■ || and inner product (■,■). Define 
A = d xx with the zero Dirichlet boundary condition. It defines a compact 
analytic semigroup S(t) , t > , on H. Denote by < Ai < A2 < • • • 
the eigenvalues of —A with the corresponding eigenf unctions {ek}^ =1 , which 
forms an orthonormal basis of H . For every a > 0, define a new norm 
\\ u \\a = || (— ^4) a//2 w||o, for those u E H such that this quantity is finite. 

Here we make the following assumptions about the mixing properties of 
the random boundary and body forcing in the random PDE system ([I])-®. 

(H g ) For every t, g{t,-) is Lipschitz continuous in u with Lipschitz con- 
stant L g and g(t,0) = 0. For every u E H, <?(-,«) is an //-valued 
stationary random process and is strongly mixing with an exponential 
rate 7 > 0. That is, 

sup sup \F(UnV) -F(U)¥(V)\ < e" 7 * , t>0, 

where < s < t < 00 , and Q\ = <r{g(r, u) : s < r <t} is the cx-algebra 
generated by {g(r, u) : s < r < t} . 

(Hf) The process f(t) is a bounded continuous differentiable process with 
1/(^)1 < Cf , for some constant Cf > , and the time derivative process 
ft(t) is a bounded stationary process with Kf t = and the mixing rate 
is exponential . That is, 

sup sup \F(U nV)- F(U)F(V) | < e~ xt , t > 0, 

where 0<s<t<oo,A>0, and J 7 * = cr{/ t (r) : s < t < t} is the 
cx-algebra generated by {ftij) : s < r <t} . 

Remark 1. A simple example of such f t is the stationary solution of the 
following linear stochastic equation 

drj = -Xr/dt + dB(t), 

where B(t) is a standard scalar Brownian motion. 
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Remark 2. Taking time derivative on the random boundary condition, we 
have 

u t = Kx + 9(t/e,u e ) , u e (x,0) = u 

u e t (0,t) = ^ft(t/e), <(/,t) = 0, 

which is a system with a random dynamical boundary condition. 

To 'homogenize' the inhomogeneous boundary condition in the system 
(fll)— ([21 , we transform the random boundary condition to the field equation 
by introducing a new random field u e = u e — y / e/(t/e)(l — f ). Then, u xx = u% x 
and the system ([I])-® becomes 

u\ = u xx + g(t/e,u<) ~ ^Jt{t/e){\ - f ), (3) 

« e (x,0) = «o->/i/(0)(l-f), (4) 
u e (0,t) = 0, u € (l,t) = 0, (5) 

which is a random system with homogeneous boundary conditions. By the 
assumption (Hy), / is bounded. Thus for every t > and x G (0, /) , 

tf-u e = \/e/(t/e)(l - f) ^ 0, e^O. (6) 

So in the following subsections, we consider u e , and derive an ensemble av- 
eraged equation to be satisfied by the limit of u € . 

2.1 Tightness 

In this section, we examine the tightness of the distribution of u e in space of 
continuous functions, C(to, T; H), for all fixed T > to > . 

In the mild or integral formulation, the equation fl3]) becomes 

u £ (t) = S(t)tf(0) + [ S(t-s)g&u e (s))ds-± [ S(t-s)f t (l)(l-f)ds. 



o Jo 



By the properties of the semigroup S(t), we have 



S(t-8)f t ($(l-£)d8 



W(t)\\o<\\V(0)\\o+j o \\g(i,u'(s))\\ods+± 

and by the assumption (H s ) , 

\\g(z,u*(s))\\ < L g \\u'(s)\\ < L g (\\u e (s)\\ + VdC f ) . (7) 

5 



Then we have, for every T > and < t < T , 

E sup ||u e (s)||o < ||'" e (0) 1 1 o + L g I sup ||u e (r)||o ds 

0<s<t Jo 0<r<s 



+ C T ,i + sup ||/ e (s)||o, (8) 

0<s<t 

where Ct,i is a positive constant depending only on L g , / and Cf , and 

P(t) :=± f S (t-s)f t (l)(l-f)ds. 
Jo 

Next we treat the singular term I e (t) . By the factorization method [2], for 
some < a < 1 , 

je^ = sir^a /" ( t - s)"' 1 S (t - S )Y £ (s) ds , 

Jo 

with 

= ^ T (s _ r)_a/t( " )5(s ■ r)(i ■ f } dr • (9) 

Then, for every T > 0, there is a positive constant CV,2 such that 
E sup \\I e (s)\\l < C T>2 [ E\\Y e (s)\\lds, 0<t<T. 

0<s<t ' Jo 

Notice that 

n\Y e (s)\\l = ± f f\ a -r)-°(8-T)-°E[M3Mz)] 
Jo Jo 

x S(s - r)(l - f)S(s - r)(l - f) drdr , 

by the assumption (H/) . For every T > , there is a positive constant Ct,3 
such that for all < t < T, 

E sup ||r( S )||o<C r , 3 . (10) 

0<s<t 

Hence, for every T > 0, applying the Gronwall inequality to (JSJ) , we obtain 
E sup ||« e (t)||o<C T (l + ||« e (0)|| ), (11) 

0<t<T 

for some constant Ct > . Furthermore, from the mild form of u e , by the 
fact that ||S'(t)«||i < 4|IM|o> we have 

< ^ll« e (0)|| + J* -^\\g(° e , U *(s))\\ ds + ll/^IK . (12) 
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We now consider the term ||/ e (s)||i . Still by the factorization method, 
||/ £ (t)||i < ^ [\t-s) a - 1 \\S(t-s)Y e (s)\\ 1 ds 



< ^ (t-s 
Jo 



,a-l 1 



't-s 



\Y e (s)\\ ds, 



where Y e (s) is defined by (Q . Then, choose a with 1/2 < a < 1, and by the 
same discussion for ffTUl) . we conclude that for every T > 

< C T , 5 , 0<t<T 

for some constant Ct,s > 0. Then for t > 0, from ([7j) and (fl2"j) . and by 
Gronwall inequality we have 



E||u ( 



< C T , t <t<T, 



(13) 



for some constant CV > . 

To show the tightness of the distributions of u e , we need a Holder estimate 
in time. For every < s < t < T , 



\\u e (t) - u e (s)\\ 
< \\(S(t)-S(s))tf(0)\\ Q + 



+ 



+ 



i 

v/i 



S{t-a)g{z,u\<j))do 
r S (t-a)f t ^)(l-f)da 

J s o 

Hsit-^-Sia-aMZu'Wda 
Jo o 

F(* -*)]/*(?)(! -f)d°- 



o 



By the estimate on ||-u e (t)||o and ([7]) , we have for some constant Ct > 
-t 



E 



S{t-a)g{%u%a))da 



< C T Vt^s~. 



o 



Moreover, by the strong continuity of the semigroup S(t), we also have 



E 



[S(t-a)-S(s-<T)]g(Zu*(a))da 
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Now consider for the singular terms. First notice that (1 — x/l) is smooth. 
We have 

S(t-a)(l-f) e L°°(0,l) . 

Therefore, 

2 



E± 



t i-t 



S(t-a)f t (i)(l-f)d* 

Js Js Jo 



t l-t 



< a 



for a positive constant C^t depending on T and I . Now by (H/), we have 



E^ 



S(t-a)f t (z)(l-Vda 



< C hT {t - s) 



Furthermore, 



[S (t-a)-S(s-a)}M^(l-f)da 

k j ° j ° 







where 



(1 - f)e fc (x) dx. 



Then still by (Hj), we have for some constant Ct > 



E^ 



< C T (t - s) 



Now we need the following lemma [6] . Suppose X\ and A^ are two Banach 
spaces. Let T>0,l<p<oo, and B be a compact operator from Ai to A2. 
That is, B maps bounded subsets of X\ to relatively compact subsets of X 2 . 

Lemma 1 ([6]). LetH be a bounded subset o/L 1 (0, T; Ai) suc/i that G = BH 
is a subset of L p (0, T; X 2 ) bounded in L r (0, T; X 2 ) with r > 1 . If 

lim + a) — u(-)\\lp(q^t;X 2 ) — uniformly for u G G , 



then G is relatively compact in L p (0,T; X 2 ) (and in C(0,T; X 2 ) if p = +oo). 
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Then by the above lemma, with X\ = Hq(0,1) and X 2 = L 2 (0,l) and B 
being the embedding from X 1 to X 2 , we have the following tightness result. 

Lemma 2. (Tightness) 

Assume that both (H g ) and (Hf) hold. For every < to < T , the distribution 
of {u t }o< e <i is tight in space C(to,T; H) . 



2.2 Ensemble averaging 



Next we use the weak convergence method [5] to pass the limit e — > . In 
this approach we construct a martingale which has the following form 



A e $(s)ds, 



for some J^-process <E>(i) defined by z\(t) and A e , which is a pseudo differ- 
ential operator to be introduced later. 

Because of the tightness of it e in space C(to, T; H) for every fixed t > , 
in order to determine the limit equation of u e in space C(to,T; H), we con- 
sider the limit of $((u e (i), </>)), for every bounded second order differentiable 
function $ : M. — > R and for every compactly supported smooth function 

<peC?(0,l). 

First, we have 

$((u e (t),^>)-$((u ,^>) 

^((u e (s),<f))(u e (s),<f xx )ds+ [ <l>X(u e (s),<p))(g{l,tf(s)),<p)ds 



i 



$'({^( S )^))(/ t (f)(l-f)^)d S . 



(14) 



To treat the singular term in (1141) . we apply a perturbation method in 
Chapter 7]. To this end, we define the following two processes 



i 



(15) 



and 



&((u e (t),p))(f t (l)(l-f),v)ds 
Then we have the following lemma. 



T 



(16) 
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Lemma 3. Assume that (Hf) holds. Then 



E\F?(t)\<CV~e, 

for some constant C > and 

E sup \Ff(t)\ ->> 0, e -»■ 0, 

0<t<T 

/or every T > . 

Proof. By the boundedness of /t and the strong mixing property, we have 

\F^(t)\<CV~e (17) 
for some constant C > . Then by the choice of $, the proof is complete. □ 

Now we apply a diffusion approximation to derive the limit of u e in the 
sense of distribution. For this we introduce the following operator 



A e $(t) =P-lim±E $(t + 5)-$(t)|7j 



(18) 



for J-'q'' 6 measurable function $(t) with sup t E|$(£)| < oo . Using Ethier and 
Kurtz's result jU Proposition 2.7.6] , we know that 



A e $(s) ds 



is a martingale with respect to J-q . Define processes Y e and Z e as follows 

Y\t) = $((^(t), <p)) - F{{t) , Z\t) = A e Y e (t) . 

A direct calculation yields 

Z\t) (19) 
= $'((n e (t), ip)) [(u e (t) : ip xx ) + (g(t/e, u e (t)), p)} 
-$"((^(t),^))((l-f)^) 2 ^/ t (^(t) 

+ &'((u e (t) : tp)) [(u e (t), <p„) + (g(l tp)] (1 - f , ^(t) 

:= Zf(f) + Z|(t) + Z|(t) . 



10 



Next we pass the limit e — > for u € (t) in space C(to, T; H) . By the conver- 
gence result of Walsh [TTJ Theorem 6.15], we only need to consider finite di- 
mensional distributions of {(u e (t), <pi), . . . , (u e (t), (p n )} for every <fii , . . . , <p n G 
C£°(0, I) . That is, we pass limit e — >• in 



Y e {t)-Y e {s)- / Z e (r)dr 



/i((u £ (ri), <^i), . . . , (u e (r n ), <p n )) } = 



E 

for every bounded continuous function h and < r\ < ■ ■ ■ < r n < T with 
any T > . Denote by u one limit point in the sense of distribution of u e 
as e — > in space C(t ,T; H) . For simplicity we assume if converges in 
distribution to u as e — > . Then by the estimates in Lemma E] we have 

Y*{t) - Y\s) -+ *««(*),*>» - H(u(s),<p)) (20) 

in distribution. 

Consider the integral term in ( I14p . First we need the following lemma 
whose proof is given in Appendix lAl 

Lemma 4. The following convergence in probability holds: 
[g(r/e, u e (r)) - g(u e (r))} dr , e -»■ . 
Then by this lemma, we have 

Z x (r)dr^{ ^(u(r),<p))[(u(r),<p xx ) + (g(u(r)),ip)]dr (21) 

J s 

in distribution as e — > . By the the estimate (fTTj) , we have 

E / |Z|(r)|dr ->•(). (22) 

J s 

Now we consider (t) . Define a bilinear operator 

(£y?,y?>=& f [\l-fMx)(l-y T My)dxdy, (23) 
Jo Jo 

where b is the variance of f t , which is constant defined as 

b:=Ef t (t)f t (t) >0. (24) 
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We again apply a perturbation method. Set 

/oo 
E[f t (l)F*(s)-±b\^]ds. (25) 

By the properties of conditional expectation and the definition of Fq , we 
have 



/oo 
E[/*(f)^(*)-§4#l ds 

/■OO /"OO 

E[/ t (f)/ t (i)-|&|jf] drds. 



Then, by the strong mixing properties of ft in the assumption (ttf) , we have 

supEF 2 e (t) = 0(e). 

t>o 

Furthermore by the same calculation as for Z e (t), we have the following 
lemma. 

Lemma 5. The following result holds: 

A^(t) = -$"«««(*), ^((l - f),^[/ t (f)F e (t) - \b] + 0(e). (26) 

Now we have the following J^-martingale 

mi ■.= mi?®,?)) - F t(t) - m) - 

$'((w e (s), if)) [{u e (s), ip xx ) + (g(u £ (s)), </?)] ds 
+ \ f\'\(tf(s)^)){Y>v^)ds + 0(e) 



By passing the limit e — > 0, the distribution of the limit u of w e solves the 
following martingale problem 



M(t) = <f>((u(t),<p))- ^({u(s),^))[(u(s)^ xx ) + {g(u(s))^)]ds 

Jo 

+ \ [ <S>"({u( S ),tp)){Zip,ip)ds, (27) 
Jo 

which is equivalent to the fact that u is the martingale solution of the fol- 
lowing stochastic PDE: 

du = [u xx + g{u)\ dt - Vb(l - f ) dB(t), (28) 
12 



where B is a usual scalar Brownian motion, and b is the variance of f t as 
defined in 

Finally, by the uniqueness of the solution to equation (128]) . we have the 
following main result on ensemble averaging under a random boundary con- 
dition. 

Theorem 1. (Ensemble averaging under a random boundary condition) 
For every to > and T > to , the solution u e , of the random PDE system (T7]|, 
converges in distribution to u in space C(t , T; H), with u solving the limit 
equation |H]). 

3 Ensemble averaging under fast oscillating 
random body forcing 

In this section, we consider the special case when the random boundary 
condition is absent. The approach to derive ensemble averaged model in the 
last section is applicable in this case. But our goal here is to further show 
that the deviation process, u e — u, can be quantified as the solution of a linear 
stochastic partial differential equation. 

We consider the following PDE with random oscillating body forcing on 
a bounded interval (0, /) 

u \ = u e xx + f(t/e,u e ), u € (x,0) = u , u e (0,t) = 0, u e (l,t) = 0. (29) 

Here we make the following assumption on the random body forcing /. 

(H) For every t , f{t,-) is continuously differentiable and Lipschitz continu- 
ous in u with Lipschitz constant Lj and f(t, 0) = . For every u G H , 
f{-,u) is an if- valued stationary random process and is strongly mixing 
with an exponential rate 7 > 0, i.e., 

sup sup \F(U n V) - F(U)¥(V) | < e" 7 ' , t > 0. 

Here s and t satisfy the condition < s < t < 00 , and = cr{/(r, u) : 
s < t < t} is the ex-algebra generated by {/(r, u) : s < r < t} . 

We introduce the notation ip (t) to quantify the mixing as follows 
<p(t) = sup sup |P(Z7 n V) - ¥(U)F(V) | . 
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By the above assumption, for any a > 

POO 

/ (p a (t)dt <oo. 
Jo 

For the random oscillating PDE (|29|) we have an averaging principle as 
above. Introduce the following averaged equation 

u t = u xx + f(u) , tt(0) = m , (30) 

where /(w) = E/(t, w) = lim^cx, i J Q T f(s/e,u) ds. Define the deviation 
process 

^(t):=^(^(t)-n(t)). (31) 
Then the following averaging principle will be established. 

Theorem 2. (Ensemble averaging under random body forcing) 
Assume that (H) holds. Then, given a T > , for every uq G if, £/ie 
solution u e (t,Uo) of 129) converges in 'probability to the solution u of (|3"0~1) in 
C(0, T; if). Moreover, the rate of convergence is \fe, that is, for any k > 
there is C£ > such that 

pj sup ||w e (t) -u(f)|| > C£y/e\ < k. (32) 

Furthermore, the deviation process z e converges in distribution in the space C(0, T 
to z, which solves the following linear stochastic PDE 

dz{t) = [z xx (t) + ~f(u(t))z(t)} dt + dW, z{0) = z{l) = 0, (33) 

where 

FW = e/1M) 

and W(t) is an H-valued Wiener process defined on a new probability space 
(Cl, J r , P) with the covariance operator 

POO 

B(u) = 2 E [(f(t,u) - /(«)) ® (/(0,«) - /(«))] dt. 
./o 

Remark 3. 27ms deviation result is similar to the averaging results for ran- 
dom PDEs in J3 Wi- 
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Proof. First by the assumption of Lipschitz property on / in (H) , and notic- 
ing that there is no singular term here, standard energy estimates yield that 
for every T > 

sup \\u^t)\\j < C T , (34) 



0<t<T 



and 



\u € (t) - u e (s)\\ < C T \t -s\, < s < t < T, 



(35) 



with some positive constant Ct • Then we have the tightness of the distribu- 
tions of u e in space C(0, T; H) for every T > . 
Notice also that u e (t) satisfies 

u e {t) = S(t)u + [ S{t- s)fXs/e,u e (s))ds, 
Jo 

and for every ip G H 

(u e (t), ip) = (S(t)u , ip) + [\s(t - s)f(s/e, u%s)), if) ds . 

Jo 

By passing the limit e — > , we can just consider the integral term in the 
above equation. By the tightness of the distributions of u e , we can follow the 
same discussion as in Appendix [A] which yields the averaged equation (|3"Dj) 
and the esitmate (1321) . 

We next consider the deviation process z e . By the definition of z e , 

? = 4. + ^Lf(*A.« e ) - /(«)] , ^ e (o) = o, 

with the zero Dirichlet boundary condition. For every a > 0, 



\\A a z* 



±_ I A « e A(t-s) [f ^ u e {s)) _f {u{s))]ds 



< 



1 / A a e A(*- S ) [/(f)U e (s)) _ /(f)U(s))]ds 



+ 



i / A«V*<*->[/(>( a )) -/(«(«))] 



:= hit) + I 2 {t). 
Notice that for < a < 1/2 , 

±\\A«e A ^[f(l,u*(s)) - f(- e Ms))}\\o < C(l + ^)L/|k e ||o, 

15 



for some constant C > . Then 

E sup h{t) < C T , 

0<t<T 

for some constant Ct > . For I 2 , by the factorization method again, we 
have 

J 3 = s»fl f (t - s ) e - 1 e A {t - s)A a Y e (s) ds, 
Jo 

where Y e is defined as 

Y*(s) = X !\s - r) e e A ^[f^,u(r)) - f(u(r))} dr. 
Jo 

Then 

E sup I 2 {t) <C T [ E|| A a Y e (s) 1 1 o ds, 

0<t<T Jo 

for some Ct > . Notice that 

1*1 f* S t* S 

\\A a Y*(s)\\l = \ (s-r) e (s-r) e A a e A ^[fC-,u(r,x))-f(u(r,x 
Jo Jo Jo 

A a e A{s ~ T) [f{^ u(r, x)) - /(u(r, x))] drdrdx . 
A standard discussion for the averaged equation yields that 

sup \\u{t)\\l<C T , 

0<t<T 

for some constant Ct > 0. Then A a e A ^ s ~ r '[f{-, u{r, x)) — f(u(r,x))] G Jq 
and A a e A ( s ~ T ' [/(-, u(r, x)) — f(u(r,x))] G J 7 )? and they are bounded real- 
valued random variables for fixed x G (0,/). Applying a mixing property jH 
Proposition 7.2.2] and choosing positive parameters a and 9 so that a + 9 < 
1/2, we have 

K\\A a Y e (s)\\l < C T , < s < T, 

and then 

E sup hit) < C T , 

0<t<T 

for some constant Ct > . So for some a > 0, 
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Furthermore, for s , t with < s < t < T , 
\\z e (t) - z^s)\\l 



2 

e 

+ 



z A{ f- r \f\r/eX(r))-f(u{r))]dr 



(I-e 



z A{s - r) [f{r/eX{r)) - Ru{r))]dr 



Then via a similar discussion as that for f )36|) . we conclude that for some 
< 7 < 1, 

E||z e (t) - z e (s)\\l < C T \t-s\\ 

which yields the tightness of the distributions of z e in C(0,T; H) . 
We decompose z £ = z\ + ^ so that 



z\ = Az\ + - /(«)] , 2^(0) = 0, 



and 



= ^ + ^[/(t/e, u e ) - /(t/e, u)\ , ^(0) = . 

For (p G C£°(0, /), we also consider the limit <&((z{(t), </?)) for every bounded 
second order differentiable function $ : R — > R in the weak convergence 
method. Notice that 



$«^(t),^»-$«0,y?» = / ^{{z\{s)^)){z\(s)^ xx )ds 

Jo 

+± r ^((^(s), ¥,»</(!, «( a )) - />(,)), ^ ds . 

Jo 

Define the following process 



/oo 
&((zt(t)^))(f(^u(t))-f(u(t))^)ds 



T 



t/e 



(37) 



A direct calculation yields that 

F*(t) :=A e ${(4(t),<p))-A'FZ(t) 
= &((zt(t),ip))(zt(t),Aip) + <f>"((zt(t),ip)) 

x i 



E[(f(t, u(t)) - /(«(*)), <p) (f(l, u(t)) - /(«(*)), | 

/■oo 
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Define two bilinear operators 

B\u, 5, t) := 2[f(l u) - /(«)] S [/(f , u) - /(«)], 

and 



:= 2 / E [(/(t, u) - /(u)) ® (/(0, «) - /»)] dt 
./o 

Then by a mixing property [H Proposition 7.2.2], we have 

E\Fi(t)\ -+0, 



E 



E 



| <5 e («(t), s, p> | ?T ds - \ (B{u{t))ip, <p 



->0, 



and 



E 



->0, 



as e — 7- . Then we also have a martingale 



./o 

-\ !\"{{z\{s)^)){B{u)^^)ds + 0{e). 



By passing the limit e — > and by the same discussion as in $21 we see that 
z\ converges in distribution to Z\, which solves 



dz! = Azi + dW, z 1 (0) = 0, 



(38) 



where W is an if-valued Wiener process defined on a new probability space 
P) with covariance operator B(u) . Furthermore, z\ converges in dis- 
tribution to Z2-, which solves 

z 2 = Az 2 + f{u)z, z 2 (0) = 0. 

Then z e converges in distribution to z with z solving f|33|) . The proof is 
complete. □ 

Remark 4. The assumption on the strong mixing property in (H) can be 
weakened as 

/■oo 

(p a (t) dt < oo, 

for some a > . In this case, we also have Theorem^ See 74 ■ f or more 
details. 
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A Proof of Lemma 14 



A similar result has been given in [10J Proposition 7]. Here we present another 
proof which gives a stronger convergence, together with the convergence rate 
in probability. 

First, under the assumption (H s ), we show that for almost all u € fl , 



[g( L £ ,q) -g{o)] dr 



(39) 



for every q e H . 
Noticing 



[g(r/e,q) - g(q)} dr G H , 



we get 



[g( z e ,q) -g(q)] dr 



= SUp j. rp 

o ipeH IMIo 



[g( z ,q) - g(q)] dr,<p 



Consider 



E/j\g^,q)-g(q)] dr,p 

E (g( L e , q) - g(q),<p) (#(f , q) - g(q), drds 



JO 



By a mixing property jU Proposition 7.2.2 ] , we have 

E/j\g(^q)-g(q)] dr,<p\ = O(e)\\cp\\ , 
which yields (|39|) . 

By the estimate in §2.11 . for every k > , there is C% > 0, which is 
independent of e, such that 



P{ \\u e (t) - u e (s) 1 1 o < C^Vt^~s) >l-K, 



(40) 



for every t > s > . Furthermore, by the tightness of the distributions of {u e } 
in space C(0, T; H) , for every n > , there is a compact set K K C C(0, T; if) 
such that 

P{u e E K K } > 1 - k . (41) 
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So we define 



{u> G f2 : events in 



(HOD and (jUJ) hold} . 



Due to the compactness of K K , for every e > 0, we only need to consider a 
finite e-net {gi ,q 2 , ■ ■ ■ , gAr} in C(0, T; H), which covers {u e }. Without loss 
of generality, we assume that qj , j = 1, 2, . . . , N , are simple functions [TO] . 

Now we consider all u E £l K . By the construction of u e and boundedness 
of / , we have for u e f2 K 



l« e (0 



»||o<C£>/r^+ViC f , 



for some constant C > . 

For every 5 > , we partition the interval [0, T] into subintervals of length 
of 5 . Then for t e [k6, (k + 1)5) , < k < [J], 



< 



+ 



^(^, M e (r))-^(r))]rfr 
b(i,« e (r))-i/(= ««(**))] dr 
[g(^ qj (kS))-g( qj (kS))] dr 
[g(u e (kS))-g(u e (r))] dr 



[g(^(k8))-g(^ qj (k5))] dr 
[g( qj (k5)) - g(u*(k5))] dr 



kd 



for some qj . Notice that, by the assumption (H 9 ) and the definition of g, g 
is also Lipschitz continuous in u with the same Lipschitz constant L g . Then 
by the assumption (H g ) and the definition of fl K , 



[g(j,u e (r)) - g(u e (r))] dr 

< T[L g C£5 + v^C + L g e + 0(y/e) + L g e + L g C^S + ^~eC] . 

Due to the arbitrary choice of 8 , e and k, and notice a similar discussion as 
that in [10], we thus complete the proof. □ 
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